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Method for Determining Position and Velocity of Targets from Signals scattered by 
the Targets ; • 

The present invention relates to a method for using signals scattered by one or more 
5 targets for determining position and velocity of each target. The invention is 

especially useful in the case where the positions in space and velocity vectors for a 
very large number of targets are to be determined. The invention is' for instance 
intended for use with the system for determining position and velocity of targets for 
which this day a patent has been applied with Hans Hellsten as inventor, and which 
10 patent application 0101661-7 is hereby incorporated by reference when regarding a 
more detailed reasoning of the system of transmitter and receiver, henceforth called 
stations. 

The invention is based on a 3-dimensional position space and a 3-dimensional 
1 5 velocity space, which together form a 6-dimensional state space with up to ~-10 20 
cells, where in principle each cell shall be assigned the value 0 (no target) or 1 
(target). The large number of cells in the state space place special requirements on 
the signal-processing method, and one possible method is dealt with in said patent 
application by Hellsten. In the present patent application an alternative method is 
20 suggested. 

Due to the large number, ~ 10 20 , of cells in the state space, known projection 
methods are slow and cumbersome. The simplest known method is where all 
conceivable target positions are formed for three stations and N targets. They will 
25 be fewer than N 3 in number, and it is easy to determine position as well as velocity 
vector for the target for every possible target position. In this way ~ N 3 candidates 
are obtained, each of which is then supported or discarded by further independent 
measurements. This gives a method of processing which requires K N 3 operations. 
If K Is large then this method can be much too slow for Interesting values of N. 

30 

A specific property of the systems for which the use of the present method Is 
intended Is. according to the above, that the information in an introductory, detecting 
phase can be regarded as binary, which can be utilised. Further the Input is from 
the beginning distributed between different stations. This means that a Ideal, 
35 distributed processing which calculates parameters of position and velocity for the 
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targets gives an advantage, as one then receives a paralleling effect on the sensor- 
near level. A comprehensive communication between stations is then necessary. 

The stumbling block is the association problem, which is to in a correct way pair 
5 together data from different radar stations. An association where it is unclear 
whether it is correct or not can be called a candidate. Gradually it becomes clear 
whether a candidate is a correct association, that is a target, or a false association, 
sometimes called a ghost The problem is then to handle the initially very large 
number of ghosts, that is to find the targets among all candidates. 

10 

The purpose of the invention is to solve this problem, which is done by the invention 
being given the features that appear from the following independent claim. Suitable 
embodiments of the invention will be evident from the remaining claims. 

15 The invention will be described in more detail in the following with reference to - 
enclosed drawings, where 

fig. 1 shows how the geometry around two stations on a plane and how 
rotation-symmetrical conceivable target positions on measurement 
20 from two stations can be transferred to points on a half-plane, 

fig. 2 shows how velocities in a 3-dimensional space can be transferred to 
non-unique velocities in said half-plane and 

fig. 3 shows how one can determine position and velocity by combining 

measurements from two pairs of stations that have a common station. 

25 

If one studies figure 1 with two stations, s 1t s 2l on the ground, which see one and the 
same target, t u in mono-static measurements, one realizes immediately that they 
place targets along a half-circle in the position space above the ground. Bi-static 
measurements between the two stations share this symmetry, so that a half-circle 

30 can be accepted or discarded in Its entirety. For each point on the half-circle, an 
estimated velocity v is determined up to a line in the velocity space. That is, two 
stations together ascribe a candidate (target or ghost) a two-surface in the state 
space. A family of such two-surfaces in the state space form a relatively compli- 
cated sub-set in the state space. It is therefore desirable with an introductory 

35 calculation where candidates can be represented by points in a reduced state space, 
rather than by areas. In the case where one studies position as well as velocity, as 
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described above, the reduced state space will consequently be 4-dimensional. 
Further some type of candidate reduction should occur in this reduced state space, 
so that the final processing in the full state space becomes simplified. 

The basic idea of the invention is to work in low-dimensional, < 6, state spaces that 
are common to several sensors. Optimally the areas in the full state space, corre- 
sponding to certain sensor registrations, will be represented by points in the lower- 
dimensional ditto. As the reduced state space is common to several sensors, one 
can under certain conditions obtain a sensor-near and local elimination of many 
ghosts by calculations which are not demanding to process. Such a reduced set of 
candidates, together with other similarly reduced sets of candidates make up input 
to a final and complete association in the full state space. This final calculation is 
then, by the relatively small size of input amounts, also surmountable in terms of 
processing intensity. 



From the symmetry in figures 1 and 2 one can see that one can allow the points on 
the arc corresponding to a candidate f, be represented by a point P(tJ in a half- 
plane P 12 on one side of the symmetry axis. In figure 2 one can see how the velocity 
v can be determined from mono-static and/or bi-static measurements from a pair of 

20 stations by Doppler measurements, except for a component in tangent to the arc 
direction. This means that the velocity lies on the line Lfs^tuv). All velocitieson 
this line will be represented by the same 2-dimensional velocity u, which is a vector 
in the half-plane P„. In reverse, it applies that if the velocity in the half-plane P 12 is 
known, one knows that the velocity lies on the line Lfs^U.u), where u is now 

25 regarded as the canonic representative of u in three dimensions. 

In contrast to the simple problem of finding the common line between measurements 
from two stations, that is the arc from the intersection of two (half-)spheres, or the 
common point between measurements from three stations, that is the point of inter- 
30 section of three (half-)spheres, the common points from two or three bi-static 
measurements between the same stations are more complicated to find, as each 
measurement generates a (harf-)ellipsoid with the stations at the focal points. It is 
more difficult to calculate common points between two ellipsoids than between 
spheres. 

35 
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The invention utilises therefore mono-static measurements between pairs of stations 
to calculate candidates. It then uses the associated bi-static measurement to 
eliminate the ghosts. In this way one does not need to calculate common points 
based on the ellipsoids, but rather looks only at the points where the mono-static 
5 measurements indicate that there are candidates and sees whether the bi-static 
measurements have one or several candidates at the distance in question. This is 
much simpler than calculating common points from bi-static measurements. In 
cases where the distance-based analysis does not give sufficient elimination of the 
number of ghosts, one can use as well the compiled Doppler information. The 
10 Doppler measurements share according to the above the symmetry of the distance 
measurements and by requiring compatible mono-static and bi-static Doppler 
measurements a further reduction in the number of ghosts is obtained. 

In the following a more detailed description is given of the method in connection with 
15 an amount of distributed I isotropic transmitters/receivers of electromagnetic of 

acoustic signals. In the concrete example we consider a number of radar stations. 
Each radar station sends out a signal that is received by the station itself and also by 
surrounding neighbours. We refer again to the above-mentioned patent application 
by Hellsten for a more detailed description of how the system can be embodied in 
20 general. 

In a first step one uses a first s, and a second s 2 transmitter/receiver of electromag- 
netic or acoustic signals and carries out mono-static measurements M 1 and M 2 from 
each transmitter/receiver, and also a bi-static measurement B 12 between the 
25 stations. First the two mono-static measurements create a number of candidates 
with 2-dimensional position and 2-dimensiona! velocity. Then these are tested 
against the result of the bi-static measurement and the candidates found in all of the 
measurements with suitable error margins are retained. 

30 In a second step one carries out corresponding calculations with one of the two 
original stations, say s,, and a third station s,. The same calculations that were 
carried out for s, and s 2 are carried out for s, and s 3 . One then sorts, for each 
calculation, conceivable targets according to the distance to the common 
transmitter/receiver, here s 1( and determines the distance bins in both calculations 

35 that simultaneously contain one or several candidates. For these candidates the 
radial Doppler velocity in the two measurements is studied. If the difference is 
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smaller than a predetemtfned value then the candidates are retained; if not then the 
are discarded. ■ . • 

Figure 3 shows how a target f, is represented by a point in the half-plane P„ and 
another in the half plane P 31 . The 3-dimensional position is given by the point of 
.ntersection of the half-circles. The true velocity vfor f, is represented by the 
velocity u u which is a vector in the plane tangent to the point P( tl ) in the half-plane 
P« and corresponding vector u 2 in the half-plane P„. Each such velocity can be 
represented according to the above on a line in the 3-dimensional space The true 
velocity v lies in the intersection between these lines Lfr.s^) and LfW,.^. 

In a third step one can carry out a third bi-static measurement B 23 between the two 
transmitters/receivers, here s 2 and s 9 , which is not the common transmitter/receiver 
according to the above. The result can be used together with the candidates that 
remam after step two above to further reduce the number of ghosts. " " 

The local treatment around three stations has through this procedure, with initially 
pa.red calculations around two stations and then including all three stations in the 
ca.culat.ons, in an elegant way requiring relatively few calculations, given a number 
of targets and possibly, as a consequence of imperfections in the measurement 
system and the application of it, some remaining ghosts. It is interesting that in 
parallel with this, in a larger system, a number of other stations can also have given 
a number of targets. The parallel handling is very favourable and an absolute pre- 
requisite according to the above for the processing to be at all possible in a iarger 
to system with many stations. 

To improve the accuracy in the calculation, and discard possible remaining ghosts 
one can of course allow further stations, which have a sufficient range to be 
included in the calculate, One can consider further mono-static as weil as further 
bi-static calculations. 

H ha, already been mentioned several times that the present invention is Intended to 
be Included in a iarger system of stations, for example that which Is presented In the 
said patent app«cation by Hans Hellsten. In this application H la stated that a 
number of station, can be placed as grid point, In an essemialty equidistant grid on 
a ,u rf aca. which limit, the aurveyed etate ,paca. e.g. a ground surface. The idea in 
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the patent application of Hellsten is that the ranges of the stations shall be such that 
each target is detected on bi-static measurement by at least six independent sets of 
transmitters and receivers. For this reason the range of the stations is 2d, if the 
distance between the stations is d p if we consider an essentially planar surface. If 
5 the surface is not essentially planar then the ranges must be adapted so that an at 
least 6-fold bi-static overlap of the surveyed state space is nevertheless achieved. 
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Claims: 

1. Method for determining position and velocity of targets from signals scattered by 
the targets, c h a ra c t e r / s e d in that one uses a first and a second station 

5 comprising transmitter/receiver (s,,s 2 ) of electromagnetic or acoustic signals and 
carries out mono-static measurements (M 1( M 2 ) from each station, and also a bi-static 
measurement between the stations (B, 2 ), and first from the two mono-static 
measurements calculates a number of target candidates with 2-dimensional position 
and 2-dimensional velocity, whereupon these are tested against the result of the bi- 

1 0 static measurement and the target candidates found in all the measurements with 
suitable error margins are retained, whereby the calculations take place in a 2- 
dimensional position space and a 2-dimensional velocity space and give the position 
and the velocity modulo rotational symmetry around a line through the two stations 
(s,,s 2 ) in that one utilises that rotational symmetry results in the ability to represent 

15 ' conceivable measured target positions in a 3-dimensional position space, which lies 
along arcs with centre on the axis of symmetry, by a point {P(t)) in a half-plane (P«) 
on one side of the axis of symmetry, which makes a 2-dimensional position space, 
and that the velocity can be represented by a velocity in a 2-dimensional velocity 
space that is a plane tangent to said point. 

20 

2. Method according to claim 1 , c h a r a c t e r / s e d in that one establishes the 
position for a target candidate in a complete 3-dimensional position space and 3- 
dimensional velocity space through corresponding calculations with either the first or 
the second station (si.Sj) and a third corresponding station (s 3 ), whereupon for each 

25 calculation one sorts target candidates according to distance to the common station 
and determines for the distances in both calculations which both contain one or 
several conceivable targets, whether the candidate shall be retained as target or 
discarded depending on the agreement of the radial velocity. 



30 



3. Method according to claim 2, c h a r a c t e r i s e d in that one carries out a 
third bi-static measurement (B u ) between the two stations (s 2l s 3 ) which are not the 
common station according to claim 2 and retains a target candidate if agreement is 
achieved, otherwise the target candidate is discarded. 



WO 02/093191 




PCT/SE02/00894 



8 



4. Method according to claim 2 or 3, characterised in that one uses 
additional stations for mono-static and/or bi-static measurements which are 
compared with the result of the earlier ones. 

5 5. System according to anyone of claims 1-4, character /sed in that said 
stations (s 1( s 2 ,s 3 ) are placed as grid points in an essentially equidistant grid on a 
surface, which limits the surveyed position space, e.g. a ground surface, with the 
distance between the grid points essentially the same size, d, and where the range 
of the signals when having an essentially planar surface are at least 2d, implying at 
10 least 6 independent bi-static configurations per grid point and that the range, in the 
case that the surface is not essentially planar, is adapted to give the same number 
of bi-static configurations as in the planar case. 
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